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Объектом исследования является система «колесо — рельс». 
Предмет — компьютерное моделирование изнашивания эле- 
ментов системы. Цель настоящей работы — изучение зависи- 
мости износа колес и рельсов от параметров кинетической 
характеристики трения. Задачи работы: проанализировать 
предлагавшиеся ранее математические функции, описываю- 
щие кинетическую характеристику трения, на наличие экс- 
тремумов, поведение (рост, убывание, монотонность, унимо- 
дальность и т. п.) в различных диапазонах изменения пара- 
метров; осуществить компьютерное моделирование движения 
колесной пары по рельсовому пути с вычислением количе- 
ственных оценок износа колес и рельсов при разных парамет- 
рах кинетической характеристики трения; установить каче- 
ственные закономерности влияния задаваемых при компью- 
терном моделировании параметров кинетической характери- 
стики на износ. При моделировании использовался программ- 
ный комплекс «Универсальный механизм». В результате по- 
лучены значения различных характеристик, связанных с изно- 
сом (мощность, суммарная и удельная работа силы трения), 
при разных сочетаниях параметров кинетической характери- 
стики трения. Результаты работы позволяют сделать следую- 
щие выводы: при компьютерном моделировании движения 
железнодорожных экипажей необходимо учитывать парамет- 
ры кинетической характеристики трения, так как в диапазонах 
значений, соответствующих реальным условиям, они могут 
существенно влиять на износ; при некоторых сочетаниях зна- 
чений вышеупомянутых параметров наблюдаются ситуации, 
когда увеличение коэффициента трения сопровождается сни- 
жением износа. Предлагается продолжить исследование в 
двух основных направлениях: установление закономерностей 
влияния на кинетическую характеристику трения и износ мно- 
гообразных реальных факторов и изучение по той же методи- 
ке влияния на износ параметров динамической (при перемен- 
ной скорости) характеристики трения. 


Ключевые слова: трение, износ, кинетическая характеристи- 
ка, скорость скольжения, коэффициент трения, удельная рабо- 
та, колесо, рельс, компьютерное моделирование, факторы. 
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Введение. Износ колес и рельсов представляет собой очень важную и актуальную проблему, от которой во 
многом зависит эффективность работы железнодорожного транспорта. Ежегодные потери в масштабе железных дорог 
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России исчисляются миллиардами рублей [1]. Основным практическим методом борьбы с износом является лубрика- 
ция, направленная на снижение трения в системе «колесо — рельс». Однако таким образом решается лишь часть про- 
блемы и далеко не всегда достигаются желаемые результаты. Большое значение имеет возможность прогнозирования 
износа и его изменений при осуществлении различных мероприятий. Такую возможность дает компьютерное модели- 
рование взаимодействия железнодорожного экипажа и пути. Одним из хорошо известных и многократно проверенных 
инструментов такого моделирования является программный комплекс «Универсальный механизм» (далее «УМ»), раз- 
работанный в Брянском государственном техническом университете коллективом под руководством профессора 
Д. Ю. Погорелова. Описание возможностей УМ представлено в [2—5]. Комплекс, в частности, позволяет моделировать 
движение экипажа на прямых и кривых участках пути, вычислять различные величины, связанные с износом колес и 
рельсов. 

Износ колес и рельсов сильно зависит от трения, на которое, в свою очередь, существенно влияют разнооб- 
разные количественные и качественные факторы. Большое значение имеет зависимость силы трения от скорости 
скольжения, поскольку качение колеса по рельсу часто сопровождается проскальзыванием. Исследование этой зави- 
симости проводилось много лет. Еще в ХХ веке на основе экспериментальных данных были получены эмпирические 
формулы, связывающие коэффициент трения и скорость скольжения. Такие исследования проводились и в ХХ веке, 
продолжаются они и в настоящее время. Различные формулы представлены в работах [6-8]. Между тем до сих пор нет 
ясности, насколько хорошо предлагаемые зависимости отражают реальные случаи, и как следует задавать различные 
их параметры в разных ситуациях с учетом многообразия действующих реальных факторов. В программном комплек- 
се УМ используется формула, которая, как будет показано ниже, по сути является частным случаем известной форму- 
лы И. В. Крагельского — С. В. Щедрова [6]. В описании комплекса УМ [5] приводятся некоторые рекомендации по 
выбору параметров в зависимости от типа локомотива, влажности, наличия лубрикации. Представляет интерес изуче- 
ние уровня чувствительности результатов компьютерного моделирования по износу к задаваемым параметрам. Кроме 
того, важно установить, каким образом изменение этих параметров в более широких пределах влияет на рассматрива- 
емые процессы. На наш взгляд, такое исследование может, с одной стороны, помочь оценить точность компьютерного 
моделирования износа колес и рельсов, а с другой стороны — расширить понимание влияния различных факторов 
зависимости трения от скорости скольжения в целом. Последняя проблема имеет в трибологии большое значение, 
поскольку, как было установлено в многочисленных теоретических и экспериментальных исследованиях, такая зави- 
симость является ключевой для возникновения фрикционных автоколебаний, которые играют значительную роль в 
различных областях техники [7]. Как показано в работах [8, 9], на эту зависимость большое влияние оказывают мно- 
гочисленные взаимосвязанные количественные и качественные факторы. Попытка оценки этого влияния методами 
планирования эксперимента была сделана в [10]. В работах [11-13] были рассмотрены особенности трибологической 
информации и выработаны специальные методы ее анализа. Ниже представлена практическая реализация одного из 
предложенных там подходов для изучения кинетической характеристики трения в задачах, связанных с износом колес 
и рельсов. Результаты сопоставлялись как с уже упомянутыми отечественными работами [9, 10], так и с зарубежными 
исследованиями [14, 15], что позволило обосновать приводимые ниже выводы. 

Анализ математических функций, выражающих кинетическую характеристику трения. Предваритель- 
ное изучение функций, выражающих кинетическую характеристику трения, выполнялось простейшими методами ма- 
тематического анализа с применением программы МАТНСАО. 

Наиболее универсальной формулой кинетической характеристики трения является формула Крагельского [7]: 

1 (»›) =(а+Ь»)е “ +а (1) 
Здесь у — скорость скольжения, а, Б, с, 4 — параметры кинетической характеристики. Физический смысл этих пара- 
метров связан сложными многофакторными зависимостями со свойствами материалов, поверхностей, конструкции 
узла трения, нагруженности и внешних условий. Качественный анализ этих зависимостей представлен в [6]. Количе- 
ственное определение а, Б, с, 4 из теоретических соображений пока не представляется возможным. В работе [7] кратко 
описана методика их экспериментального определения. Подробно различные подходы проанализированы в [8]. 
В программном комплексе УМ заложена более простая по сравнению с (1) формула: 


8()=^((1-А)е +4) @) 
Очевидно, что формула (2) является частным случаем формулы (1), т. е. функция } (у) переходит в & (у), если: 
4= А; с=В; а= (1-4) $ 50. 


Качественный физический смысл параметров А и В наглядно проявляется при поочередном приравнивании их 
к некоторым характерным значениям. Если А = 0, то получаем: 


8()=Ле”. (а) 
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Формула (2а) выражает монотонно убывающую функцию, и параметр В определяет скорость убывания. Ве- 
личина В существенно зависит от температуры и давления. 
Если А =1, то 


8()=Л, (25) 

т. е. кинетическая характеристика вырождается в константу — коэффициент трения покоя. 
Если В = 0, то при любом А приходим к случаю (25), т. е. при равенстве параметра В нулю функция 8(у) также 
вырождается в константу, являющуюся коэффициентом трения покоя. Если В находится в интервале (—0;0), то 


функция 8(») является возрастающей. Если параметр В находится в интервале (0; +0 ), то функция является убываю- 
щей. 
Найдем первую и вторую производные функции (у) по у. После преобразований получим: 
' —Ву 
8'(>)=(А-ПВе”, (3) 
8®) = (1-А) В?е №. (4) 
Согласно [7], в трибологических задачах (в частности, для взаимодействия колеса и рельса), как правило, 
О<А<1иВ>0. Анализ формул (2)-—(4) показывает, что такая ситуация соответствует монотонно убывающей функ- 
ции (у), и ее график является вогнутым. 
Рассмотрим теперь функцию }„р(А, В), т. е. ту же кинетическую характеристику трения, задаваемую функцией 


8(^) по формуле (2), но при неизменной скорости и переменных параметрах А и В. С целью упрощения значение ко- 
эффициента трения }„„ = $, получаемое по формуле (2), делится на константу / (обозначим частное как }р2: 








а ое +А. (5) 
Л 
При В = соп5Е 
Вир В В В 
——=(е "-Ае "+А)'= АЦ-е ^^”), 6 
94 ( ) =А( ) (6) 
д 
ти > | при В >ОиА > 0. 
ОА 
При А = соп5Е 
Гир —_ (р-Ву\! —Ву\! — _ —Ву —Ву —Ву — 
=(е ")'-А(е ')'=-Ве '+АВе ”” =Ве "(А-\), (7) 
ОВ 
д 
Чт приВ >ОиА > 1. 
ОВ 
Ор 
При А < 1, Е <0, т. е. с ростом В ]„„з уменьшается. 


Методы, инструменты и результаты компьютерного моделирования. В качестве инструмента для прове- 
дения компьютерного моделирования использовался программный комплекс УМ — в частности, его модуль ОМ 
Госо, позволяющий моделировать движение железнодорожных экипажей. Был выбран компонент «колесная пара с 6 
степенями свободы». Рассматривалось движение на прямых и кривых участках, причем брались правые кривые и со- 
ответственно вычислялся износ левых колес в колесной паре (по умолчанию в ОМ Г.осо рассматриваются правые кри- 
вые, поэтому износ обычно определяется на левых (набегающих) колесах). В данной работе использовалась упрощен- 
ная расчетная схема: одиночная колесная пара, катящаяся по рельсовому пути. В результате моделирования определя- 
лись количественные характеристики, непосредственно связанные с износом: мощность силы трения (М»»„), работа 
силы трения (А»„„,), удельная работа силы трения (5„.„). Мощность силы трения при двухточечном контакте опреде- 
ляется как скалярное произведение силы крипа (трения) на гребне и скорости колес относительно рельса в точке кон- 
такта и рассчитывается по формуле: 


ЛИ неа = —Р 2х — Ру, (Вт), (8) 
1 
Арана | М ее (Дж), (9) 
0 
Аа» 
нех т ж/м), (10) 


где 5 — пройденный путь. 
Некоторые задававшиеся параметры моделирования представлены на рис. 1. 
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#=Ю[(1-А)ехр(-В*/5) +А] 
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Рис. 1. Параметры моделирования изнашивания элементов системы «колесо — рельс» 


Во вкладке программы «Колесо/Рельс — Трение» (рис. 1) задаются коэффициенты А и В, входящие в иссле- 
дуемую зависимость силы трения от скорости скольжения. Здесь же приведена формула (2), используемая в програм- 
ме. 
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Отложить переменную по абсциссе 



































Отложить "время" по абсциссе 


Рис. 2. Задание параметров зависимости трения от скорости 


Результаты моделирования — значения 5 „и,, А лег И М „ог, вЫЧИСЛенные по формулам (8)—(10) по найденным 
для данных моментов времени силовым и кинематическим переменным, можно перенести в М5 Ехсе! (фрагмент 
представлен на рис. 3). Столбец «Абсцисса» содержит модельное время в секундах, в других столбцах содержатся 
значения параметров в указанный момент времени. 
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АРВ СЕНО Е | 
1 Абсцисса 5$\Меаг1_11 Абсцисса А\М/еаг1_11 Абсцисса М\/еаг1_11 
236) 4.68 5.2779846 4.68 493.847534 4.68 267.419434 
237 4.7 5.3125162 4.7 499.2005 4.7 267.856842 
238 4.72 5.346848 4.72 504.562195 4.72 268.287109 
239 4.74 5.380981 4.74 509.932373 4.74 268.710693 
240 4.76 5.4149156 4.76 515.310974 4.76 269.128174 
241 4.78 5.4486542 4.78 520.697876 4.78 269.539032 
242] 4.8 5.4821968 4.8 526.092896 4.8 269.943848 
243 4.82 5.5155454 4.82 531.495972 4.82 270.342865 
244 4.84 5.5487008 4.84 536.906982 4.84 270.735748 
245 4.86 5.5816641 4.86 542.325806 4.86 271.123138 


Рис. 3. Фрагмент результатов в окне М$ Ехсе] 


Компьютерное моделирование производилось многократно с различными значениями параметров АиВ. На 
рис. 4 представлены графики, построенные на основании полученных результатов. По оси абсцисс указаны значения 
параметра В в диапазоне от —1 до 1, а по оси ординат — значения, связанные с износом исследуемых величин при ука- 
занных выбранных значениях параметра А. 
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Фактор износа, мощность 
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Рис. 4. Графики характеристик износа при некоторых значениях А 


На рис. 5 представлены графики изменения тех же связанных с износом величин при некоторых различных 
значениях В. По оси абсцисс здесь отложены значения параметра А. 
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Рис. 5. Графики характеристик износа при некоторых значениях В 


Шаг и диапазоны изменения по параметрам А и В выбирались, исходя из намерений изучить чувствитель- 
ность результатов к их изменению. Надо особо отметить, что ряд расчетов выполнялся с отрицательными значениями 
В в соответствии с возрастающей кинетической характеристикой трения. В реальном взаимодействии колеса и рельса 
такие случаи встречаются относительно редко, но они возможны при некоторых условиях смазывания и, кроме того, 
представляют теоретический интерес. 

В табл. 1 представлены данные, полученные в результате анализа рис. 4, 5 и характеризующие разброс иссле- 
дуемых параметров износа. 

Таблица 1 


Сводные показатели износа 





тт тах | рост, % 
5 иеаг 3,38 7,85 132,45 

| Ам | 338,90 | 787,90 | 132,46 | 
Му 124,40 | 485,90 | 290,62 


























Как видно из табл. 1, разброс (разница между наименьшим и наибольшим значением) значительный по всем 
трем характеристикам, особенно чувствительной оказалась величина М „еа»- 

В табл. 2 содержатся значения коэффициентов А и В, приведенные непосредственно в программном комплек- 
се УМ в качестве предлагаемых при компьютерном моделировании для некоторых типов локомотивов и при различ- 
ном состоянии рельсов (влажные или сухие). 

Таблица 2 


Рекомендуемые в УМ значения коэффициентов А и В 












































Наименование А В 
Локомотив, сухой рельс 0,4 0,6 
Локомотив, влажный рельс 0,4 0,2, 
Локомотив $045Х, сухой рельс 0,44 0,6 
Локомотив ОВ127, сухой рельс 0,38 0,7 
Локомотив $252, сухой рельс 0,36 0,55 
Локомотив 3ВВ460, влажный рельс 0,5 0,16 
Локомотив 12Х, влажный рельс 0,4 0,4 
Локомотив $045Х, влажный рельс 0,38 0,18 
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Анализ разброса результатов, собранных в табл. 3, относящихся к предлагаемым в УМ значениям А и В, пока- 
зывает, что в этой области разница между наибольшим и наименьшим значениями характеристик износа не превыша- 
ет 22 %. Надо отметить, что сама область является весьма узкой, и многие факторы, влияющие на кинетическую ха- 
рактеристику трения (температура окружающей среды, наличие и тип смазки, загрязнения и др.) здесь не учтены. Эти 
вопросы требуют дополнительного всестороннего исследования. 

Таблица 3 


Сводные показатели по предустановленным в программе параметрам 





тт тах | рост, % 
О НЕВЕ 4,14 4,80 16 
Анеат 415,07 | 481,94 16 
Мьеаь 159,40 | 194,18 22 





























По результатам компьютерного моделирования можно сделать также следующие выводы: 
—Щ при увеличении значения параметра В увеличиваются значения исследуемых характеристик износа уе Анеаь» 
М иеаг» 
—Щ при увеличении значения коэффициента А уменьшаются значения исследуемых характеристик износа ая А неа» 
М хеаь» кроме случая при В= 0, когда эти характеристики не изменяются. 

Отметим, что сопоставление последнего вывода со значениями найденных выше первых частных 

Вир о Вир 

ОА ОВ 
сопровождаться уменьшением работы силы трения и износа. 


Ои 








производных (при В >Ои0<А<1 <0), говорит о том, что увеличение коэффициента трения может 


Выводы. Результаты проведенного исследования позволяют сделать следующие выводы. 

1. В области значений параметров кинетической характеристики трения, соответствующих реальным ситуа- 
циям взаимодействия колеса и рельса, влияние их изменений на износ в общем случае является существенным, т. е. 
при компьютерном моделировании взаимодействия колеса и рельса надо учитывать зависимость коэффициента тре- 
ния от скорости, если в задачу входит оценка износа. 

2. Как показывает компьютерное моделирование движения колесной пары по рельсовому пути, при падаю- 
щей кинетической характеристике трения уменьшение коэффициента трения может сопровождаться увеличением ха- 
рактеристик износа. 

3. В трибодинамических системах значение коэффициента трения является функцией мгновенной скорости 
скольжения, которая изменяется во времени. Поэтому для компьютерного моделирования движения железнодорож- 
ных экипажей желательно использовать не кинетическую, а динамическую характеристику трения, которая может от 
нее существенно отличаться. Это требует проведения дополнительных теоретических и экспериментальных исследо- 
ваний. 
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